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ABSTRACT 

An all sky map of the apparent temperature and optical depth of thermal dust emission is constructed using the Planck-HFI (350 jjm to 2 mm) 
and IRAS (100//m) data. The optical depth maps are correlated with tracers of the atomic (Hi) and molecular gas traced by CO. The correlation 
with the column density of observed gas is Unear in the lowest column density regions at high Galactic latitudes. At high A^h, the correlation is 
consistent with that of the lowest A^h, for a given choice of the CO-to-H2 conversion factor. In the intermediate A^h range, a departure from linearity 
vN is observed, with the dust optical depth in excess of the correlation. This excess emission is attributed to thermal emission by dust associated with a 
dark gas phase, undetected in the available H i and CO surveys. The 2D spatial distribution of the dark gas in the solar neighbourhood (|bii| > 10°) 
is shown to extend around known molecular regions traced by CO. 

The average dust emissivity in the Hi phase in the solar neighbourhood is found to be Td/NJ^^ = 5.2 x 10~^^ cm^ at 857 GHz. It follows roughly a 
power law distribution with a spectral index j3 - l.S all the way down to 3 mm, although the SED flattens slightly in the millimetre. Taking into 
• • account the spectral shape of the dust optical depth, the emissivity is consistent with previous values derived from FIRAS measurements at high 
^ latitudes within 10%. The threshold for the existence of the dark gas is found at NjJ^ = (8.0+0.58)xlO^^ Hcm~^ (Ay = 0.4 mag). Assuming the same 
high frequency emissivity for the dust in the atomic and the molecular phases leads to an average Xco = (2.54 + 0.13) X 10^^ H2cm~^/(Kkm s~^). 
The mass of dark gas is found to be 28% of the atomic gas and 118% of the CO emitting gas in the solar neighbourhood. The Galactic latitude 
distribution shows that its mass fraction is relatively constant down to a few degrees from the Galactic plane. 
A possible explanation for the dark gas lies in a dark molecular phase, where H2 survives photodissociation but CO does not. The observed 
transition for the onset of this phase in the solar neighbourhood (Ay = 0.4 mag) appears consistent with recent theoretical predictions. It is also 
possible that up to half of the dark gas could be in atomic form, due to optical depth effects in the Hi measurements. 

Key words. ISM: general, dust, extinction, clouds - Galaxies: ISM - Infrared: ISM - Submillimeter: ISM 



1. Introduction 

The matter that forms stars, that is left over after star formation, 
or that has never experienced star formation comprises the in- 

* Corresponding author; email: Jean-Philippe.Bernard@cesr.fr. 



terstellar medium (ISM). The life-cycle and the duration of the 
various observable phases remains largely unknown, because the 
nature of the diffuse interstellar medium is difficult to discern, 
owing to its low temperatures and large angular scales. 
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The distribution of diffuse interstellar gas, by which we mean 
gas not in gravitationally-bound structures and not in the im- 
mediate vicinity of active star-formation regions, has primarily 
been assessed using the 21 -cm hyperfine line of atomic hydro- 
gen. That line is easily excited by collisions and is optically 
thin for gas with temperature Tk > 50 K and velocity disper- 
sion > lOkms"^ as long as the column density is less than 
9 X 10^^ cm"^ (Kulkarni & Heiles 1988). Such conditions are 
typical of the diffuse ISM pervaded by the interstellar radiation 
field (ISRF), because photoelectric heating from grain surfaces 
keeps the gas warm (T > 50 K), and observed velocity disper- 
sions (presumably due to turbulence) are typically > lOkms"^ 
Based on the observed dust extinction per unit column density, 
N(HI)/Av = 1.9 X 10^1 cm-2 mag-^ (Bohlin et al. 1978), the 
upper column density for optically thin 21 -cm lines corresponds 
to visible extinctions Ay < 4.7. Thus the 21 -cm line is expected 
to trace diffuse, warm atomic gas accutately throughout the dif- 
fuse ISM, except for lines of sight that are visibly opaque or are 
particularly cold. 

Molecular gas is typically traced by the 2.6-mm 
^^CO(/=1^0) rotational line in emission, which, like the 
21 -cm Hi line, can be easily excited because it involves energy 
levels that can be obtained by collisions. The CO emission line, 
however, is commonly optically thick, due to its high radiative 
transition rate. In the limit where the lines are optically thick, the 
primary information determining the amount of molecular gas 
in the beam is the line width. If the material is gravitationally 
bound, then the virial mass is measured and CO can be used 
as a tracer of molecular mass. It is common astronomical 
practice to consider the velocity-integrated CO line intensity 
as measuring the molecular column density, with the implicit 
assumption that the material is virialized and the mass of the 
virialized structures is being measured. In the diffuse ISM, 
these conditions typically do not apply. On a physical scale of R 
(measured in parsecs), interstellar material is only virialized if 
its column density N > 5.2 x lO^^^V^R"^ cm"^ where is the 
velocity dispersion (measured in kms"^). Thus the diffuse ISM 
is typically gravitationally unbound, invalidating the usage of 
CO as a virial tracer of the molecular gas mass, except in very 
compact regions or in regions that are visibly opaque. Although 
CO can emit in gas with low density, the critical density required 
for collisional equilibrium is of order 10^ cm"^, which further 
complicates the usage of CO as a tracer. This again is not typical 
of the diffuse ISM. 

To measure the amount and distribution of the molecular 
ISM, as well as the cold atomic ISM, other tracers of the in- 
terstellar gas are required. At least three tracers have been used 
in the past. These are UV absorption in Werner bands of H2, 
infrared emission from dust, and y-ray emission from pion pro- 
ductiondue to cosmic-rays colliding with interstellar nucleons. 
The UV absorption is exceptionally sensitive to even very low 
H2 column densities of 10^^ cm"^. Using Copernicus (Savage 
et al. 1977) and FUSE data, atomic and molecular gas could 
be measured simultaneously on the sightlines to UV-bright stars 
and some galaxies. A survey at high Galactic latitudes with 
FUSE showed that the molecular fraction of the ISM, f(H2) = 
2N(H2)/[2N(H2) + N(HI)] < 10"^ for lines of sight with total 
column density less than 10^^ cm"^, but there is a tremendous 
dispersion from 10""^ to 10"^ for higher-column density lines of 
sight (Wakker 2006). Since UV-bright sources are preferentially 
found towards the lowest-extinction sightlines, an accurate aver- 
age f(H2) is extremely difficult to determine from the stellar ab- 
sorption measurements. Along lines of sight toward AGNs be- 
hind diffuse interstellar clouds, Gillmon & Shull (2006) found 



molecular hydrogen fractions of 1-30% indicating significant 
molecular content even for low-density clouds. 

The dust column density has been used as a total gas column 
density tracer, with the assumption that gas and dust are well 
mixed. The possibility that dust traces the column density bet- 
ter than H i and CO was recognized soon after the first all- sky 
infrared survey by IRAS , which for the first time revealed the 
distribution of dust on angular scales down to 5'. Molecular gas 
without CO was inferred from comparing IRAS 100 jim surface 
brightness to surveys of the 21 -cm and 2.6-mm lines of Hi and 
CO on 9' or degree scale de Vries et al. (1987); Heiles et al. 
(1988); Blitz et al. (1990). At y scale using Arecibo, the cloud 
G236-h39 was found to have significant infrared emission un- 
accounted for by 21 -cm or 2.6-mm lines, with a large portion 
of the cloud being possibly H2 with CO emission below detec- 
tion threshold (Reach et al. 1994). Meyerdierks & Heithausen 
(1996) also detected IR emission surrounding the Polaris flare in 
excess of what was expectated from the H i and CO emission, 
which they attributed to diffuse molecular gas. The all sky far- 
infrared observations by COBE -DIRBE (Hauser et al. 1998) 
made it possible to survey the molecular gas not traced by H i 
or CO at the 1° scale (Reach et al. 1998). This revealed numer- 
ous "infrared excess" clouds, many of which were confirmed 
as molecular after detection of faint CO with NANTEN (Onishi 
et al. 2001). Finally, there are also indications of more dust emis- 
sion than seen in nearby external galaxies such as the Large 
Magellanic Cloud (Bernard et al. 2008; Roman-Duval et al. 
2010) and the Small Magellanic Cloud (Leroy et al. 2007). This 
suggests that large fractions of the gas masses of these galaxies 
are not detected using standard gas tracers. 

The y-rays from the interstellar medium provide an indepen- 
dent tracer of the total nucleon density. As was the case with the 
dust column density, the y-ray inferred nucleon column density 
appears to show an extra component of the ISM not associated 
with detected 21 -cm or 2.6-mm emission; this extra emission 
was referred to as "dark gas" (e.g. Grenier et al. 2005; Abdo et al. 
2010), a term we will adopt in this paper to describe interstellar 
material beyond what is traced by H i and CO emission. Grenier 
et al. (2005) inferred dark gas column densities of order 50% of 
the total column density toward locations with little or beyound 
detection threshold CO emission, and general consistency be- 
tween infrared and y-ray methods of detection. Recent observa- 
tions using FERMI have significantly advanced this method, al- 
lowing y-ray emission to be traced even by the high-latitude dif- 
fuse ISM. In the Cepheus, Cassiopeia, and Polaris Flare clouds, 
the correlated excess of dust and y rays yields dark gas masses 
that range from 40 % to 60 % of the CO-bright molecular mass 
(Abdo et al. 2010). 

Theoretical work predicts a dark molecular gas layer in re- 
gions where the balance between photodissociation and molec- 
ular formation allows H2 to form in significant quantity while 
the gas-phase C remains atomic or ionized (Wolfire et al. 2010; 
Glover et al. 2010). In this paper we describe new observa- 
tions made with Planck ^ (Planck Collaboration 201 la) that trace 
the distribution of submillimeter emission at 350 yum and longer 
wavelengths. In combination with observations up to lOOyum 
wavelength by IRAS and COBE -DIRBE , we are uniquely able 



^ Planck (hiip://www.essL.mi/Planck) is a project of the European 
Space Agency (ESA) with instruments provided by two scientific con- 
sortia funded by ESA member states (in particular the lead countries: 
France and Italy) with contributions from NASA (USA), and telescope 
reflectors provided in a collaboration between ESA and a scientific con- 
sortium led and funded by Denmark. 
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Table 1. Characteristics of the data used in this study 



LJaVa 


^ref 


Vref 


U 


O' II 


O' abs 




[yum] 


[GHz] 


[arcmin] 


[MJy/sr] 


[%] 


IRAS 


100.0 


2998 


4.30 


0.06*^* 


13.6* 


HFI 


349.8 


857 


3.67 


0.12^ 


7% 


HFI 


550.1 


545 


3.80 


0.12^ 


7% 


HFI 


849.3 


353 


4.43 


0.08^ 


<2% 


HFI 


1381.5 


217 


4.68 


0.08^ 


<2% 


HFI 


2096.4 


143 


7.04 


0.08^ 


< 2% 


HFI 


2997.9 


100 


9.37 


0.07^ 


<2% 


Data 


line 


^ref 


e 


cr 


abs 






[//m] 


[arcmin] 


[Kkms-i] 


[%] 


LAB 


Hi 


21 cm 


36.0 


1.70^ 


10.0 


DHT 


i^CO 


2.6 mm 


8.8 


1.20 


24.0 


Dame 


i^CO 


2.6 mm 


8.4 


0.6 


24.0 


NANTEN 


i^CO 


2.6 mm 


2.6 


1.20 


10.0 



^ Assumed to be for the average IRAS coverage, ctu computed by rescal- 
ing this value to actual coverage. * From Miville-Deschenes & Lagache 
(2005). Icr average value in one beam scaled from Planck HFI Core 
Team (201 lb). We actually use internal variance maps for ctu ^ Icr aver- 
age value. We actually use a map of the uncertainties (see Sect. 2.2.1). 



to trace the distribution of interstellar dust with temperatures 
down to ~ 10 K. The surface brightness sensitivity of Planck, 
in particular on angular scales of 5' to 7°, is unprecedented. 
Because we can measure the dust optical depth more accurately 
by including the Planck data, we can now reassess the relation- 
ship between dust and gas, and relate it to previous infrared and 
independent UV and y-ray results, and compare it to theoretical 
explanations to determine just how important the dark gas is for 
the evolution of the interstellar medium. 



2, Observations 

2.1. Planck data 

The Planck first mission results are presented in Planck 
Collaboration (2011a) and the in-flight performances of the two 
focal plane instruments HFI (High Frequency Instrument) and 
LFI (Low Frequency Instrument) are given in Planck HFI Core 
Team (2011a) and IMennella et al. (2011) respectively. The data 
processing and calibration of the HFI and LFI data used here 
is described in Planck HFI Core Team (2011b) and Planck 
Collaboration (2011b) respectively. 

Here we use only the HFI (DR2 release) data, the process- 
ing and calibration of which are described in Planck HFI Core 
Team (2011b). In this data the CMB component was identified 
and subtracted through a Needlet Internal Linear Combination 
(NILC) (Planck HFI Core Team 201 lb). 

We use the internal variance on intensity (cTjj) estimated dur- 
ing the Planck data processing and provided with the Planck- 
HFI data, which we assume represents the white noise on the 
intensity. Note that this variance is inhomogeneous over the sky, 
owing to the Planck scanning strategy (Planck Collaboration 
2011a), with lower values in the Planck deep fields near the 
ecliptic poles. We have checked that, within a small factor (< 2), 
the data variance above is consistent with "Jack-Knife" maps 
obtained from diff'erencing the two halves of the Planck rings. 
We also use the absolute uncertainties due to calibration uncer- 
tainties given in Planck HFI Core Team (2011b) for HFI and 



summarized in Table 1 . We note that, for a large scale analysis 
such as carried out here, variances contribute to a small fraction 
of the final uncertainty resulting from combining data over large 
sky regions, so that most of the final uncertainty is due to abso- 
lute uncertainties. 

2.2. Ancillary data 

2.2.1. HI data 

In order to trace the atomic medium, we use the LAB 
(Leiden/ Argentine/Bonn) survey which contains the final data 
release of observations of the Hi 21 -cm emission line over 
the entire sky (Kalberla et al. 2005). This survey merged the 
Leiden/Dwingeloo Survey (Hartmann & Burton 1997) of the 
sky north at ^ > -30° with the lAR (Instituto Argentino de 
Radioastronomia) Survey (Amal et al. 2000; Bajaja et al. 2005) 
of the Southern sky at ^ < -25°. The angular resolution and the 
velocity resolution of the survey are ~ 0.6° and ~ 1.3 km s"^ 
The LSR velocity range -450 < Vlsr < 400 km s"^ is fully cov- 
ered by the survey with 891 channels with a velocity separation 
of AVch = 1.03 kms-^ 

The data were corrected for stray radiation at the Institute for 
Radioastronomy of the University of Bonn. The rms brightness- 
temperature noise of the merged database is slightly lower in the 
southern sky than in the northern sky, ranging over 0.07-0.09 K. 
Residual uncertainties in the profile wings, due to defects in the 
correction for stray radiation, are for most of the data below a 
level of 20 to 40 mK. We integrated the LAB data in the ve- 
locity range -400 < Vlsr < 400 km s"^ to produce an all sky 
map of the H i integrated intensity (Wri), which was finally pro- 
jected into the HEALPix pixelisation scheme using the method 
described in Sect. 2.3.1. 

We estimate the noise level of the Whi map as 
ATrmsAVch ^|Nch where Nch(= 777) is the number of channels 
used for the integration, and AT^is is the rms noise of the in- 
dividual spectra measured in the emission-free velocity range 
mainly in -400 < Vlsr < 350 km s"^. The resulting noise of 
the Whi map is mostly less than ~ 2.5 Kkms"^ all over the sky 
with an average value of ~ 1.7 Kkms"\ except for some limited 
positions showing somewhat larger noise (~ 10 Kkms"^). 

2.2.2. CO data 

In order to trace the spatial distribution of the CO emission, we 
use a combination of 3 large scale surveys in the ^^CO(/=1^0) 
line. 

In the Galactic plane, we use the Dame et al. (2001) survey 
obtained with the CfA telescope in the north and the CfA-Chile 
telescope in the south, referred to here as DHT (Dame, Hartmann 
& Thaddeus). These data have an angular resolution of 8.4' ±0.2' 
and 8.8' ±0.2' respectively. The velocity coverage and the veloc- 
ity resolution for these data vary from region to region on the sky, 
depending on the individual observations composing the survey. 
The most heavily used spectrometer is the 500 kHz filter bank 
providing a velocity coverage and resolution of 332 km s"^ and 
1.3 kms"\ respectively. Another 250 kHz filter bank providing 
the 166 kms"^coverage and 0.65 km s"^ resolution was also fre- 
quently used . The rms noises of these data are suppressed down 
to 0.1-0.3 K (for details, see their Table 1). The data cubes have 
been transformed into the velocity-integrated intensity of the line 
(Wco) by integrating the velocity range where the CO emis- 
sion is significantly detected using the moment method proposed 
by Dame (2011). The noise level of the Wco map is typically 
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-3.0 ^^^^^^^^^ ^^^^^H 30.0 K>^km/s 



Fig. 1. Map of the ^^CO(/=1^0) integrated intensity used in this paper combining the Dame et al. (2001) and high latitude survey 
and the NANTEN survey. The data shown cover 62.8% of the sky. The map is shown at a common resolution of all the sub-surveys 
of 8.8'. Many small clouds at high latitude are not visible in this rendering of the data. 



-1.2 Kkms ^ , but it varies by a factor of a few depending on the 
integration range used. 

We also use the unpublished high latitude survey obtained 
using the CfA telescope (Dame et al. 2010, private communica- 
tion). This survey is still on-going and covers the northern sky 
up to latitudes as high as |bn| = 70° which greatly increases 
the overall sky coverage. The noise level of the CO spectra are 
suppressed to ~ 0.18 K for the 0.65 km s" Velocity resolution, 
and the total CO intensity was derived by integrating typically 
10-20 velocity channels, which results in a noise level of 0.4- 
0.6Kkms-i. 

Finally, we combined the above survey with the NANTEN 
^^CO(/=1^0) survey obtained from Chile. This survey comple- 
ments some of the intermediate Galactic latitudes not covered by 
the Dame et al. (2001) maps with an angular resolution of 2.6'. 
Most of the survey along the Galactic plane has a velocity cov- 
erage of ~ 650 km s"^ with a wide band spectrometer, but a part 
of the survey has a coverage of ~ 100 km s"^ with a narrow band 
spectrometer. The noise level achieved was 0.4-0.5 K at a ve- 
locity resolution of 0.65 kms"^ The CO spectra were sampled 
with the 2' grid in the Galactic centre, and with the 4' and 8' 
grid along the Galactic plane in the latitude range |b| < 5 °and 
|b| > 5°, respectively. The integrated intensity maps were ob- 
tained by integrating over the whole velocity range, excluding 
regions of the spectra where no emission is observed. The re- 
sulting rms noise in the velocity-integrated intensity map varies 
depending on the width of the emission. This survey along the 
Galactic plane is still not published in full, but parts of the survey 
have been analyzed (e.g. Fukui et al. (1999); Matsunaga et al. 
(2001); Mizuno & Fukui (2004)). A large amount of the sky at 



intermediate Galactic latitude toward the nearby clouds is also 
covered with a higher velocity resolution of ~ 0.1 km s"^ with a 
narrow band spectrometer with a < 100 km s"^ band (e.g. Onishi 
et al. (1999); Kawamura et al. (1999); Mizuno et al. (2001)). 
The velocity coverage, the grid spacing, and the noise level for 
these data vary, depending on the characteristics of the individ- 
ual clouds observed, but the quality of the data is high enough to 
trace the total CO intensity of the individual clouds. 

The three surveys were repixelised into the HEALPix pix- 
elisation scheme (Gorski et al. 2005) with the appropriate pixel 
size to ensure Shannon sampling of the beam (Nside=2048 for 
the NANTEN2 survey and Nside=1024 for the CfA surveys) us- 
ing the procedure described in Sect. 2.3.1. 

Each survey was smoothed to a common resolution of 
8. 8' through convolution with a Gaussian with kernel size ad- 
justed to go from the original resolution of each survey to a 
goal resolution of 8.8', using the smoothing capabilities of the 
HEALPix software. We checked the consistency of the differ- 
ent surveys in the common region observed with NANTEN and 
CfA. 

We found a reasonably good correlation between the two but 
a slope indicating that the NANTEN survey yields 24% larger 
intensities than the CfA values. The origin of this discrepancy 
is currently unknown. We should note that the absolute inten- 
sity scale in CO observations is not highly accurate as noted 
often in the previous CO papers. Since the CfA survey covers 
most of the regions used in this paper and has been widely used 
for calibrating the H2 mass calibration factor Xco, in particular 
by several generations of gamma ray satellites, we assumed the 
CfA photometry when merging the data, and therefore rescaled 
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the NANTEN data down by 24% before merging. Note that this 
an arbitrary choice. The implications on our results will be dis- 
cussed in Sec. 6.1. 

The 3 surveys were then combined into a single map. In do- 
ing so, data from different surveys falling into the same pixel 
were averaged using cr^ as a weight. The resulting combined 
map was then smoothed to the resolution appropriate to this 
study. The resulting CO integrated intensity map is shown in 
Fig.l. 

2.2.3. IRdata 

We use the IRIS (Improved Reprocessing of the IRAS Survey) 
IRAS lOOyum data (Miville-Deschenes & Lagache 2005, see) in 
order to constrain the dust temperature. The data, provided in 
the original format of individual tiles spread over the entire sky 
were combined into the HEALPix pixelisation using the method 
described in Sect. 2.3.1 at a HEALPix resolution (Nside = 2048 
corresponding to a pixel size of 1.7'). The IRAS coverage maps 
were also processed in the same way. We assume the noise prop- 
erties given in Miville-Deschenes & Lagache (2005) and given 
in Table 1. The noise level of Q.06MJy sr~^ at 100 yum was as- 
sumed to represent the average data noise level and was appro- 
priately multiplied by the coverage map to lead to the pixel vari- 
ance of the data. 

2.3. Additional Data processing 

2.3.1 . Common angular resolution and pixelisation 

The individual maps are then combined into HEALPix using the 
intersection surface as a weight. This procedure was shown to 
preserve photometry accuracy. 

The ancillary data described in Sect. 2.2 were brought to the 
HEALPix pixelisation, using a method where the surface of the 
intersection between each HEALPix pixel with each FITS pixel 
of the survey data is computed and used as a weight to regrid 
the data. The HEALPix resolution was chosen so as to match 
the Shannon sampling of the original data at resolution 0, with 
a HEALPix resolution set so that the pixel size is < ^/2.4. The 
ancillary data and the description of their processing will be pre- 
sented in Paradis & et. al. (2011). 

All ancillary data were then smoothed to an appropriate res- 
olution by convolution with a Gaussian smoothing function with 
appropriate FWHM using the smoothing HEALPix function, 
and were brought to a pixel size matching the Shannon sampling 
of the final resolution. 

2.3.2. Background levels 

Computing the apparent temperature and optical depth of ther- 
mal dust over the whole sky requires accurate subtraction of 
any off'set (/^) in the intensity data, either of instrumental or 
astrophysical origin. Although both the IRIS and the Planck- 
HFI data used in this study have been carefully treated with 
respect to residual off'sets during calibration against the FIRAS 
data, the data still contains extended sources of emission un- 
related to the Galactic emission, such as the Cosmic InfraRed 
Background (CIB) signal (Miville-Deschenes et al. 2002; Planck 
Collaboration 201 Id) or zodical light which could aff'ect the de- 
termination of the dust temperature and optical depth at low sur- 
face brightness. 

In order to estimate the above data off'sets, we first compute 
the correlation between IR and Hi emission in a reference re- 



Table 2. Thermal Dust emissivity derived from the correlation 
with Hi emission in the reference region with |bn| > 20° and 
N^^ < 1.2 X Hcm-2 ((I^/Nh)'^^. Off'sets derived from an 
empty region with N^^ < 2.0 x 10^^ Hcm"^, assuming the same 
emissivity (ly/Nn)^^^. The uncertainties are given at the Icr level. 
The corresponding data are plotted in Fig. 2. 



y 




offset 


(GHz) 


[MJy/sr/lO^o Hcm-^] 


[MJy/sr] 


IRAS 


2998 


(6.95+0.94)xl0-i 


(7.36+0.03)xl0-i 


DIRBE : 


2998 


(6.60+0.01)xl0-i 


(7.72+0.22)xl0-i 


2141 


1.16+0.01 


1.41+0.16 


1249 


(8.85+0.04)xl0-i 


(8.46+0.90)xl0-i 


Planck-HFI : 


857 


(5.43+0.38)xl0-i 


(2.57+0.05)xl0-i 


545 


(1.82+0.13)xl0-i 


(1.83+0.05)xl0-i 


353 


(4.84+0.10)xl0-2 


(9.50+0.24)xl0-2 


217 


(1.14+0.03)xl0-2 


(3.45+0.1 l)xl0-2 


143 


(2.92+0.07)xl0-3 


(1.01+0.06)xl0-2 


100 


(1.13+0.04)xl0-3 


(3.34+0.60)xl0-3 


Planck-LFI : 


70.3 


(8.74+2.66)xl0-5 


(-6.40+0.93) 10-3 


44.1 


(1.14+0.16)xl0-4 


(-7.02+0.45) 10-3 


28.5 


(2.11+0.14)xl0-4 


(-6.95+0.15)10-3 


WMAP: 


93.7 


(-1.05+0.57)xl0-4 


(-1.65+1.02)xl0-3 


61.2 


(-1.14+0.27)xl0-4 


(-4.90+3.63)xl0-4 


41.1 


(3.52+1. 13)xl0-5 


(9.76+1.28)xl0-4 


32.9 


(1.39+0.08)xl0-4 


(4.10+0.97)xl0-4 


23.1 


(2.69+0.04)xl0-4 


(8.49+0.44)xl0-4 



gion such that |bn| > 20° and < 1.2 x 10^^ Hcm"^. This was 
done using the IDL regress routine and iterative removal of out- 
liers. The derived dust emissivities ((ly/Nn)^^^) are given in Table 
2. The uncertainties given are those derived from the correla- 
tion using the data variance as the data uncertainty. The derived 
emissivities are in agreement with the ensemble average of the 
values found for the local H i velocities in Planck Collaboration 
(201 li) (see their Table 2) for individual smaller regions at high 
Galactic latitude, within the uncertainties quoted in Table 2. Note 
that these emissivities are used only to derive the off'sets in this 
study. 

We then select all sky pixels with minimum H i column den- 
sity defined as < 2.0 x 10^^ Hcm"^ and compute the average 
H I column density in this region to be N^°^^ = 1.75x10^^ Hcm"^. 
The off'sets are then computed assuming that the dust emissivity 
in this region is the same as in the reference region, ie, 

/^ = /^*-(/JArHrfxN^* (1) 

where tf^^ is the average brightness in the hole region at fre- 
quency y. 

The off'set values derived from the above procedure are given 
in Table 2 and were subtracted from the maps used in the rest of 
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Fig. 2. Upper panel: Thermal Dust emissivity (Iv/Nh) from 
Table 2. The dot curve showns a grey body at Td = 17.5 ^ and 
13= 1.8 normalized at 857 GHz, for comparison. The various 
colours are for different instruments: IRAS (Yellow), DIRBE 
(light blue), Planck-HFI (red), WMAP (dark blue) and Planck- 
LFI (green). Lower panel: Offsets from Table 2. The error bars 
are plotted at ±3(T. 



this analysis. The offset uncertainties also listed in Table 2 were 
derived from the emissivity uncertainties propagated to the offset 
values through Eq. 1. When subtracting the above offsets from 
the IRAS and Planck intensity maps, the data variances were 
combined with the offset uncertainties in order to reflect uncer- 
tainty on the offset determination. Note that, for consistency and 
future use. Table 2 also lists emissivities and offset values for 
FIR-mm datasets not used in this study. Note also that these off- 
sets for Planck data are not meant to replace the official values 
provided with the data, since they suppress any large scale emis- 
sion not correlated with H i, whatever their origin. 

3. Dust temperature and emissivity 

3.1. Temperature determination 

As shown in previous studies (e.g. Reach et al. 1995; Finkbeiner 
et al. 1999; Paradis et al. 2009; Planck Collaboration 2011e,i), 
the dust emissivity spectrum in our Galaxy cannot be represented 
by a single dust emissivity index /S over the full FIR-submm 
domain. The data available indicate that JS is usually steeper in 
the FIR and flatter in the submm band, with a transition around 
500 yum. As dust temperature is best derived around the emission 
peak, we limit the range of frequencies used in the determina- 



tion to the FIR, which limits the impact of a potential change of 
13 with frequency. 

In addition, the dust temperature derived will depend on the 
assumption made about JS, since these two parameters are some- 
what degenerate in;^^ space. In order to minimize the above ef- 
fect, we derived dust temperature maps using a fixed value of 
the dust emissivity index JS. The selected JS value was derived by 
fitting each pixel of the maps with a modified black body of the 
form ly oc v^By(TD) in the above spectral range (method referred 
to as "free yS"). This leads to a median value of Td = 17.7 K and 
JS = 1.8 in the region at |bn| > 10°. Note that the JS value is con- 
sistent with that derived from the combination of the FIRAS and 
Planck-HFI data at low column density in Planck Collaboration 
(20111). Inspection of the corresponding Td and JS maps indeed 
showed spurious values of both parameters, caused by their cor- 
relation and the presence of noise in the data, in particular in low 
brightness regions of the maps. 

We then performed fits to the FIR emission using the fixed 
= 1.8 value derived above (method referred to as "fixed yS"). 
In the determination of Td, we used the IRIS lOOyum map and 
the two highest HFI frequencies at 857 and 545 GHz. Although 
the median reduced is slightly higher than for the "free yS" 
method, the temperature maps show many fewer spurious val- 
ues, in particular in low brightness regions. This results in a 
sharper distribution of the temperature histogram. Since we later 
use the temperature maps to investigate the spectral distribution 
of the dust optical depth and the dust temperature is a source of 
uncertainty, we adopt the "fixed yS" method maps in the follow- 
ing. The corresponding temperature and uncertainty maps are 
shown in Fig. 3. 

Temperature maps were derived at the common resolution of 
those three channels as well as at the resolution of lower inten- 
sity data. The model was used to compute emission in each pho- 
tometric channels of the instruments used here (IRAS , Planck- 
HFI ), taking into account the colour corrections using the ac- 
tual transmission profiles for each instrument and following the 
adopted flux convention. In the interest of computing efficiency, 
the predictions of a given model were tabulated for a large set 
of parameters (Tj^^jS). For each map pixel, the;^^ was computed 
for each entry of the table and the shape of the distribution 
around the minimum value was used to derive the uncertainty on 
the free parameters. This included the effect of the data variance 
(Tjj and the absolute uncertainties. 

3.2. Angular distribution of dust temperature 

The all-sky map of the thermal dust temperature computed as 
described in Sec. 3.1 for /S = 1.8 is shown in Fig. 5. The elon- 
gated regions with missing values in the map correspond to the 
IRAS gaps, where the temperature cannot be determined from 
the Planck-HFI data alone. The distribution of the temperature 
clearly reflects the large scale distribution of the radiation field 
intensity. 

Along the Galactic plane, a large gradient can be seen from 
the outer Galactic regions, with Td - 14 - 15 K to the in- 
ner Galactic regions around the Galactic center regions with 
Td - 19^. This asymmetry was already seen at lower angular 
resolution in the DIRBE (Sodroski et al. 1994) and the FIRAS 
(Reach et al. 1995) data. The asymmetry is probably due to the 
presence of more massive stars in the inner Milky Way regions, 
in particular in the molecular ring. The presence of warmer dust 
in the inner Galaxy is actually clearly highlighted by the radial 
distribution of the dust temperature derived from Galactic inver- 
sion of IR data (e.g. Sodroski et al. 1994; Paladini et al. 2007; 
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Fig. 3. Upper panel: All sky map of the dust temperature in K. The temperature is derived from modeling the IRIS lOOyum and 
the Planck-HFI emission at 857 and 545 GHz. Lower panel: All sky map of the dust temperature uncertainty in %. The maps are 
shown in Galactic coordinates with the Galactic centre at the centre of the image. Grey regions correspond to missing IRAS data. 



Planck Collaboration 201 If). The origin of the large scale region but we note that it corresponds to a region of enhanced X-ray 
near (ln,bn)=(340°,-10°) with Td - 20K is currently unclear, emission in the Rosat All-Sky Survey (RASS). 
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Fig. 4. Details of the dust temperature (left column), dust optical depth at 857 GHz (central column) and dark gas column density 
(right column) for the Chamaeleon (first line), Aquila-Ophiuchus flare (second line), Polaris flare (third line) and Taurus (fourth 
line). The temperature and optical depth maps are shown in log scale with a colour scale ranging from 15 K to 20 K and 1 x 10"^ 
to 3 X 10"^ respectively. The dark gas column density derived from the optical depth at 857 GHz (see Sec. 4) and is shown in linear 
scale with a colour scale ranging from -3 to 7 x 10^^ Hcm"^ . The contours show the ^^CO(/=1^0) integrated intensity at 2, 10 and 
20 Kkms"^ The double line shows the limit of the CO surveys. 



It may therefore correspond to warm dust associated with hot 
gas pervading the local bubble around the Sun, or a pocket of hot 
gas in Loop I. Similar large regions with enhanced dust temper- 
ature, such as near (lii,bii)=(340°,-30°) or (ln,bii)=(315°,+30°) 
may have a similar origin. Loop I (ln,bn)=(30°,+45°) is seen as 
a slightly warmer than average structure at Td - 19 K. Running 
parallel to it is the Aquila-Ophiuchus flare (lii,bii)=(30°,-h20°) 
with apparent Td - 14 K extending to latitudes as high as 



60°. The Cepheus and Polaris Flare (Iii,bii)=(100-120°,-Fl0- 
-h20°) (see Planck Collaboration 20111, for a detailed study) is 
also clearly visible as a lower temperature arch extending up to 
bn=30° into the North Celestial Pole loop and containing a col- 
lection of even colder condensations (Td - 12 - 13 K). 

On small angular scales, which are accessible over the whole 
sky only with the combination of the IRAS and Planck-HFI 
data at 5', the map shows a variety of structures that can all 
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be identified with local heating by known single stars or Hii 
regions for warmer spots and with molecular clouds for colder 
regions. Figure 4 illustrates the high resolution spatial distribu- 
tion of dust temperature and dust optical depth around some 
of these regions. Warmer regions include the tangent direc- 
tions to the spiral Galactic arms in Cygnus (ln,bn)=(80°,0°) and 
Carina (ln,bn)=(280°,0°), hosts to many OB associations, and 
many H ii regions along the plane. At higher Galactic latitude, 
dust heated by individual hot stars such in the Ophiuchi region 
(lii,bii)=(340°,-h20°) with individual stars cr - Sco, v - Sco, 
p - OpK ^ - OpK in Orion (ln,bii)=(210°,-20°) with the 
Trapezium stars or in Perseus-Taurus (ln,bn)=(160°,-20°) with 
the CaHfornia Nebula (NGC1499) can clearly be identified. Note 
the Spica HII region at (lii,bii)=(300°,-h50°) where dust temper- 
atures are Td - 20 ^ due to heating by UV photons from the 
nearby (80 pc) early-type, giant (BlIII) star a Vir. 

At intermediate and high latitudes, nearby molecular clouds 
generally stand out as cold dust environments with Td - 
13 K. The most noticeable ones are Taurus (ln,bn)=(160°,-20°) 
(see Planck Collaboration 2011j, for a detailed study), RCrA 
(ljj,bii)=(0°,-25°), Chamaeleon (ln,bii)=(300°,-20°) and Orion 
(libbn)=(200°,-20°). Numerous cold small scale condensations 
can readily be found when inspecting the temperature map, 
which mostly correspond to cold cores similar to those dis- 
covered at higher resolution in the Herschel data (e.g. Andre 
et al. 2010; Konyves et al. 2010; Molinari et al. 2010; Juvela 
et al. 2010) and in the Planck Cold-Core catalog (see Planck 
Collaboration 2011g,h). 

Individual nearby Galaxies are also readily identified, in par- 
ticular the Large (ln,bii)=(279°,-34°) and the Small Magellanic 
Cloud (lii,bii)=(30r ,-44°) (see Planck Collaboration 201 Ic, for 
a detailed study), as well as M31 and M33. 

Near the Galactic poles, the temperature determination be- 
comes noisy at the 5 'resolution due to the low signal levels. 

3.3. Optical depth determination 

Maps of the thermal dust optical depth (td(/1)) are derived using: 

Iv(^) 



td(^) = 



Bv(Td)' 



(2) 



where By is the Planck function and ly{X) is the intensity map 
at frequency v. We used resolution-matched maps of Td and 
Iy(/1) and derived Td(/1) maps at the various resolutions of the 
data used here. The maps of the uncertainty on td(/1) (Atd) are 
computed as: 



Atd(v) = Td 



6^y 



(Td) 



B2(Td) J 



.1/2 



(3) 



4. Dust/Gas correlation 

We model the dust opacity (tm) as 



tm(/1) 



(\ref 



[Nhi + 2XcoWco], 



(4) 



where is the reference dust emissivity measured in low 

A^H regions and Xco = /Wco is the traditional H2/CO con- 
version factor. It is implicitly assumed that the dust opacity per 
unit gas column density is the same in the atomic and molecular 
gas. If this is not the case, this will directly impact our derived 



Xco since only the product of Xco by the dust emissivity in the 
CO phase can be derived here. The fit to derive the free 

parameters of the model is performed only in the portion of the 
sky covered by all surveys (infrared, H i, and CO) and where ei- 
ther (1) the extinction is less than a threshold Ay^, or (2) the CO 
is detected with Wco > 1 Kkms"^ Criterion (1) selects the low- 
column density regions that are entirely atomic and suff'er very 
small H I optical depth eff'ects, so that the dust in this region will 
be associated with the Hi emission at 21 -cm. Criterion (2) se- 
lects regions where the CO is significantly detected and the dust 
is associated with both the Hi and the ^^CO emission lines. We 
fit for the following three free parameters '•{^) , Xco and Ay^ . 

The threshold Ay^ measures the extinction (or equivalently the 
column density) where the correlation between the dust optical 
depth and the Hi column density becomes non-linear. 

The correlation between the optical depth for various photo- 
metric channels and the total gas column density (N{^^ = Nhi + 
2XcoWco) is shown in Fig. 6. The correlations were computed 
in the region of the sky where the CO data is available (about 
63% of the sky) and at Galactic latitudes larger than bn> 10°. 
The Td and Wco maps used were smoothed to the common res- 
olution of the Hi map (0.6°). For these plots, we used a fixed 
value of Xco = 2.3 x 10^^ H2cm"^/(Kkm s"^). The colours show 
the density of points in N{^^ and td bins. The dots show the N{^^ 
binned average correlation. The larger scatter of these points at 
high N{^^ comes from the limited number of points in the corre- 
sponding bins. The red line shows the tm model values derived 

from the fit (slope=(^) ) to the low N{^^ part of the data. 

It can be seen that the correlation is linear at low N{^^ values 
and then departs from linear at N{^^ ^ 8.0 x 10^^ Hcm"^ (A^^ ^ 
0.4 mag). Above Nj^^ ^ 5 x lO^^ Hcm'^ (A^^^^^ ^ 2.5 mag), 
where N{^^ becomes dominated by the CO contribution, the dust 
optical depth again is consistent with the observed correlation at 
low N{^^ for this given choice of the Xco value. Between these 
two limits, the dust optical depth is in excess of the linear cor- 
relation. The same trend is observed in all photometric channels 
shown, with a similar value for the threshold. It is also observed 
in the HFI bands at lower frequencies, but the increasing noise 
at low Nj^^ prevents an accurate determination of the fit parame- 
ters. 

The best fit parameters for , Xco and Ay^ are given 

in Table 3. They were derived separately for each frequency. 
The uncertainty was derived from the analysis of the fitted 

(\ref 
values decrease with increas- 
ing wavelength, as expected for dust emission. The resulting dust 
optical depth SED is shown in Fig. 7. The dust optical depth 
in low column density regions is compatible with p = 1.8 at 
high frequencies. The best fit p value between the IRAS 100 /im 
and the HFI 857GHz_is actually found to be /5 = 1.75. The 
SED then flattens slightly at intermediate frequencies with a 
slope of p = 1.57 around A = 500 yum then steepens again to 
/3 = 1.75 above 1mm. The Xco values derived from the fit 
are constant within the error bars, which increase with wave- 
length. The average value, computed using a weight propor- 
tional to the inverse variance is given in Table 3 and is found 
to be Xco = 2.54 ± 0.13 x 10^^ H2cm-V(Kkms-i). Similarly, 
the 0.4 parameter does not significantly change over the whole 
frequency range and the weighted average value is found to be 
A^^ = 0.4 ± 0.029 mag. 

The excess column density is defined using the diff'erence be- 
tween the best fit and the observed dust opacity per unit column 
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Fig. 5. Maps of the dust optical depths on a log scale, in the IRAS lOOyum (first row left) and Planck-HFI bands at 857 (first 
row right), 545 (second row left), 353 (second row right), 217 (third row left), 143 (third row right) and 100 GHz (fourth row). 
All maps are shown in Galactic coordinates with the Galactic centre at the centre of the image. The missing data in all images 
correspond to the IRAS gaps. The upper and lower bounds of the colour scale are set to Tmin = 5 x 10"^ x (/1/lOOyum)"^-^ and 
Tmax = 10"^ X (/I/ lOOyum)"^"^ respectively. 



density using. 



Nh/ 



ref 



(5) 



The map is used to derive the total excess mass (M^) assum- 
ing a fiducial distance to the gas responsible for the excess. 

We also computed the atomic and molecular total gas masses 
over the same region of the sky, assuming the same distance. In 
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Fig. 6. Correlation plots between the dust optical depth at IRAS 100 yum (upper left), HFI 857 GHz (upper right), 545 GHz (lower 
left) and 353 GHz (lower right) and the total gas column density N^^^ in the solar neighbourhood (|bn| > 10°). The color scale 
represents the density of sky pixels on a log scale. The blue dots show a N^^^ -binned average representation of the correlation. The 

red line shows the best linear correlation derived at low N^^^ values ('^ = (^) * N^^^ + cste). The vertical lines show the positions 
corresponding to Ay = 0.37 mag and Ay = 2.5 mag. These figures are shown for a single Xqq = 2.3 x 10^^ H2cm"^/(Kkm s"^). 



the region covered by the CO survey, the H i to CO mass ratio de- 
rived for Xco = 2.54 x 10^^ H2cm-V(Kkm s"^) is Mhi/Mco=4. 

Using the average and Xco values above, the ratio of 

the dark gas mass to the atomic gas mass (M^/M™) and to the 
molecular gas mass (M^/M^^) are given in Table 3. On average, 
at high Galactic latitudes, the dark gas masses are of the order of 
28%% ± 3%% of the atomic gas mass and ^ 118%% ± 12%% 
of the molecular mass. 

5. Dark-gas spatial distribution 

The spatial distribution of the dark gas as derived from td com- 
puted from the HFI 857 GHz channel is shown in Fig. 8. It 
is shown in the region where the CO data are available and 
above Galactic latitudes of |bnl > 5°. Regions where Wco > 
1 Kkm/s have also been excluded. The map clearly shows that 
the dark gas is distributed mainly around the best known molec- 
ular clouds such as Taurus, the Cepheus and Polaris flares, 
Chamaeleon and Orion. The strongest excess region is in the 
Aquila-Ophiuchus flare, which was already evident in Grenier 
et al. (2005). 

Significant dark gas is also apparent at high latitudes, south 
of the Galactic plane in the anticenter and around known translu- 



cent molecular clouds, such as MBM53 (lii= 90°, bii= -30°). 
As with all the molecular clouds, the spatial distribution of the 
dark gas closely follows that of the Gould-Belt (Perrot & Grenier 
2003) and indicates that most of the dark gas in the solar neigh- 
bourhood belongs to this dynamical structure. 

6. Discussion 

6.1. Dust emissivity in the atomic neutral gas 

In the solar neighbourhood, Boulanger et al. (1996) measured 
an emissivity value in the diffuse medium of 10"^^ cm^/H at 
250 yum assuming a spectral index /3 = 2 which seemed con- 
sistent with their data. The optical depth of dust derived in 
our study in the low N{^^ regions at |bn| > 10° is shown in 
Fig. 7. The Figure also shows the reference value by Boulanger 
et al. (1996) which is in good agreement with the values de- 
rived here, interpolated at 250 yum (in fact 10% above when us- 
ing = 1.8 and 6% above when using j3 = 1.75). Our study 
does not allow us to measure the emissivity in the molecular 
gas, since we are only sensitive to the product of this emissivity 
with the Xco factor. However, we note that our derived average 
Xco = 2.54 X 10^^H2cm"^/(Kkms"^) is significantly higher 
than previously derived values. Even if we account for the pos- 
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observational results, we must put them into the same units. We 
define the dark gas fraction as the fraction of the total gas column 
density that is dark, 



A (/xrr 



Fig. 7. Dust optical depth derived from this study using the 
IRAS and Planck-HFI frequencies. The square symbol shows 
the emissivity at 250 jim derived by Boulanger et al. (1996). The 
dash and dash-dot lines show a power law emissivity with A~^-^ 
and A~^-^^ respectively, normalized to the data at lOOyum. The 
error bars shown are ±lcr. 



sible uncertainty in the calibration of the ^^CO(/=1^0) emis- 
sion (24%) discussed in Sec. 2.2.2, increasing the CO emission 
by the corresponding factor would only lower our Xco estimate 
to Xco = 2.2 H2cm"^/(Kkms"^). In comparison, a value of 
(1.8 ± 0.3) X 1020H2cm-V(Kkms-i) was found at |bn| > 5° 
from the comparison of the Hi, CO, and IRAS lOOyum maps 
(Dame et al. 2001). Similarly, values derived from y-ray FERMI 
data can be as low as Xco = 0.87 x 10^^ H2cm"^/(Kkms"^) in 
Cepheus, Cassiopea and Polaris (Abdo et al. 2010). This could 
be evidence that the dust emissivity in the high-latitude molecu- 
lar material could be larger than in the atomic phase by a factor 
^ 3. Such an increase in the dust emissivity in molecular regions 
has been inferred in previous studies (e.g. Bernard et al. 1999; 
Stepnik et al. 2003) and was attributed to dust aggregation. 

6.2. Dark Molecular gas 

The nature of 'dark molecular gas' has recently been investi- 
gated theoretically by Wolfire et al. (2010), who specifically ad- 
dress the HI/H2 and C/C^ transition at the edges of molecular 
clouds. The nominal cloud modeled in their study is relatively 
large, with total column density 1.5 x 10^^ cm"^, so the applica- 
bility of the results to the more translucent conditions of high- 
Galactic-latitude clouds is not guaranteed. The envelope of the 
cloud has an H i column density of 1 .9 x 10^^ cm"^ which is more 
typical of the entire column density measured at high latitudes. 
Wolfire et al. (2010) define foe as the fraction of molecular gas 
that is dark, i.e. not detected by CO. In the nominal model, the 
chemical and photodissociation balance yields a total H2 col- 
umn density of 7.0 x 10^^ cm"^, while the 'dark' H2 in the tran- 
sition region where CO is dissociated has a column density of 
1.9 X 10^^ cm"^. The fraction of the total gas column density that 
is molecular, 

^^"^^ - 2N(H2) + N(HI) 

is 93% in the nominal model, which suggestss that the line of 
sight through such a cloud passes through material which is al- 
most entirely molecular. To compare the theoretical model to our 
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2A^(Hf^^) 
2A^(H2) + N(HI) 



^ /(H2)fDG. 



(7) 



For the nominal Wolfire et al. (2010) model, /dg =0.29 so we 
can infer /dark=0-27. The smaller clouds in Figure 11 of their 
paper have larger /dg, but /(H2) is also probably smaller (not 
given in the paper) so we cannot yet definitively match the model 
and observations. These model calculations are in general agree- 
ment with our observational results, in that a significant fraction 
of the molecular gas can be in CO-dissociated 'dark' layers. 

If we assume that all dark molecular gas in the solar neigh- 
bourhood is evenly distributed to the observed CO clouds, the 
average foo measured is in the range foo = 1.06 - 1.22. This is 
more than three times larger than predicted by the Wolfire et al. 
(2010) mass fraction. This may indicate that molecular clouds 
less massive than the ones assumed in the model actually have a 
dark gas mass fraction higher by a factor of about three. This 
would contradict their conclusion that the dark mass fraction 
does not depend on the total cloud mass. 

The location of the Hi-to-H2 transition measured here 
(Ay^ ^ 0.4 mag) is comparable, although slightly higher than 
that predicted in the Wolfire et al. (2010) model (A^^ ^ 
0.2 mag). Again, this diff'erence may indicate variations with the 
cloud size used, since UV shadowing by the cloud itself is ex- 
pected to be less efficient for smaller clouds, leading to a transi- 
tion deeper into the cloud. 

6.3. Other possible origins 

The observed departure from linearity between td and the ob- 
servable gas column density could also in principle be caused by 
variations of the dust/gas ratio (D/G). However, such variations 
with amplitude of 30% in the solar neighbourhood and a sys- 
tematic trend for a higher D/G ratio in denser regions would be 
difficult to explain over such a small volume and in the presence 
of widespread enrichment by star formation. However, the fact 
that the dark gas is also seen in the y-ray with comparable ampli- 
tudes is a strong indication that it originates from the gas phase. 
The dark gas column-densities inferred from the y-ray observa- 
tions are also consistent with a standard D/G ratio (Grenier et al. 
2005). 

The observed excess optical depth could also in principle be 
due to variations of the dust emissivity in the FIR-Submm. We 
expect such variations to occur if dust is in the form of aggre- 
gates with higher emissivity (e.g. Stepnik et al. 2003) in the dark 
gas region. We note however that such modifications of the op- 
tical properties mainly affect the FIR-submm emissivity and are 
not expected to modify significantly the absorption properties in 
the Visible. Therefore, detecting a similar departure from linear- 
ity between large-scale extinction maps and the observable gas 
would allow us to exclude this possibility. 

Sky directions where no CO is detected at the sensitivity of 
the CO survey used (0.3-1.2 Kkms"^) may actually host signifi- 
cant CO emission, which could be responsible for the excess dust 
optical depth observed. Evidences for such diffuse and weakly 
emitting CO gas have been reported. For instance, in their study 
of the large-scale molecular emission of the Taurus complex. 
Goldsmith et al. (2008) have found that half the mass of the com- 
plex is in regions of low column density Nh < 2 x 10^^ cm"^, 
seen below Wco - 1 Kkms"^ However, Barriault et al. (2010) 
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Fig. 8. Map of the excess column density derived from the 857 GHz data. The map is shown in Galactic coordinates with the 
Galactic centre at the centre of the image. The grey regions correspond to those where no IRAS data are available, regions with 
intense CO emission (Wco > 1 Kkms"^) and the Galactic plane (|bnl < 5°). 



Table 3. Derived parameters for the dark gas, computed in the region with available CO data and |bn| > 10°. 



Frequency 
[GHz] 


[10-25 cm^] 


[H2cm-V(Kkms-0] 


[mag] 


Mi/M^i 




2998 


4.66+0.65 


(2.60+0.18)xl020 


(4.05+0.39)xlO-i 


(2.91+0.38)xlO-i 


1.27+0.16 


857 


(5.25+0.49)xl0-i 


(2.52+0.29)xl020 


(3.92+0.64)xl0-i 


(2.73+0.66)xl0-i 


1.22+0.30 


545 


(2.53+0.29)xl0-i 


(2.52+0.36)xl020 


(3.96+0.79)xl0-i 


(2.77+0.95)xl0-i 


1.24+0.43 


353 


(1.18+0.17)xlO-i 


(2.3 l+0.47)x 10^0 


(4.07+1.28)xl0-i 


(2.26+0.67)xl0-i 


1.10+0.33 


217 


(6.03+1. 19)xl0-2 


(2.58+0.67)xl020 


(4.60+2.39)xl0-i 


(1.71+4.60)xl0-i 


(7.46+20.09)xl0-i 


143 


(2.98+1. 15)xlO-2 


(1.83+1.00)xl020 


(4.62+4.81)xl0-i 






100 


(2.08+0.48)xl0-2 


(4.05+ 1.93)x 10^0 


(6.69+5. 15)xlO-i 






Average 




(2.54+0.13)xl020 


(4.03+0.29)xl0-i 


(2.78+0.28)xl0-i 


1.18+0.12 



reported a poor spatial correlation between emission by dif- 
fuse CO and regions of FIR excess in two high Galactic lati- 
tude regions in the Polaris Flare and Ursa IVLajor. The difficulty 
at finding the CO emission associated to dark gas is that the 
edges of molecular coulds tend to be highly structured spatially, 
which could explain why many attempts have been unsuccess- 
ful (see for instance Falgarone et al. 1991). In our case, it is 
possible to obtain an upper limit to the contribution of weak 
CO emission below the survey detection threshold, by assum- 
ing that pixels with undetected CO emission actually emit with 
Wco = 0-5 Kkms"^ This is the detection limit of the survey we 
use at \b\ > 10° so this should be considered an upper limit to the 
contribution of undetected diff'use CO emission. In that case, the 
dark gas mass is reduced by a factor lower than 20%. This indi- 
cates that, although diff'use weak CO emission could contribute 



a fraction of the observed excess emission, it cannot produce the 
bulk of it. 

Finally, we recognize that the optically thin approximation 
used here for the Hi emission may not fully account for the 
whole atomic gas present, even at high latitude. H i emission is 
subject to self absorption and Nh can be underestimated from 
applying too high a spin temperature (Ts) while deriving column 
densities. Ts is likely to vary from place to place depending on 
the relative abundance of CNIM clumps (with thermodynamical 
temperatures of 20-100 K) and WNIM clouds (at several thou- 
sand K) in the telescope beam. The eff'ective spin temperature 
of 250-400 K to be applied to correct for this blending andto 
retrieve the total column density from the H i spectra does not 
vary much in the Galaxy (Dickey et al. 2003, 2009). It indi- 
cates that most of the Hi mass is in the warm phase and that 
the relative abundance of cold and warm Hi is a robust frac- 
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0.0 t ^ \ ^ ^ ^ \ ^ ^ ^ \ ^ ^ ^ \ ^ ^ ^ ^ ^ ^ ^ \ ^ 

2 4 6 8 10 

b^i^ (degrees) 

Fig. 9. Fractional mass of the dark gas with respect to the neutral 
gas mass as a function of the lower bn value used in the analysis. 
The solid curve is computed under the assumption of optically 
thin H I, the dashed curve is for N[J^ computed using Ts = 80 K. 
Error bars are Icr. 



tion across the Galaxy (outside of the inner molecular ring). The 
correlation between the FERMI y-ray maps and the H i column 
densities derived for different spin temperatures also support an 
average (uniform) effective spin temperature > 250^ on and 
off the plane (Ackermann et al. 2010). In order to test these ef- 
fects, we performed the analysis described in this paper using a 
very low choice for the H i spin temperature. We adopted a value 
of Ts = 80 K when the observed H i peak temperature is below 
80KandTs = 1.1 xT 

peak when above. Under this hypothesis, we 
obtained dark gas fractions which are about half of those given 
in Table 3 under the optically thin approximation. We consider 
this to indicate that significantly less than half of the detected 
dark gas could be dense, cold atomic gas. We further note that, 
under the optically thin Hi hypothesis, the dark gas fraction ap- 
pears very constant with Galactic latitude down to |bn| - 3° (see 
Sec. 6.4), while it varies more strongly using Tg = 80 K. This 
does not support the interpretation that the bulk of the dust ex- 
cess results from underestimated Hi column densities. 

6.4. Dark-Gas variations with latitude 

We investigate the distribution of the dark gas as a function of 
Galactic latitude. This is important, since the dark gas template 
produced here for the solar neighbourhood is also used in direc- 
tions toward the plane for Galactic inversion purpose in Planck 
Collaboration (201 If). We performed the calculations described 
in Sec. 4 for various values of the Galactic latitude lower cutoff 
(bmin) in the range b^in < |bn| < 90° with bmin varying from 0° 
to 10°. For each value, we used the best fit parameters derived 
from bmin = 10° and given in Table 3. 

Figure 9 shows the evolution of the dark gas mass fraction 
with respect to the atomic gas mass as a function of bmin- It can 
be seen that the ratio changes only mildly (increases by a factor 
1.12 from hmin = 10° to bmin = 2°) as we approach the Galactic 
plane. This indicates that a fairly constant fraction of the dark 
gas derived from the solar neighbourhood can be applied to the 
rest of the Galaxy. 

Figure 9 also shows the same quantity computed using the 
H I column density derived using Ts = 80 K. It can be seen that. 



in that case, the dark gas fraction is predicted to decrease by a 
factor 2.12 from bmin = 10° to bmin = 2°. This is caused by 
the much larger inferred H i masses toward the plane under this 
hypothesis. We consider it unlikely that the dark gas fraction 
varies by such a large factor from the solar neighbourhood to the 
Galactic plane, and consider it more likely that the correction 
applied to Nh by using a spin temperature as low as Tg = 80 K 
actually strongly overestimates the H i opacity, and therefore the 
fraction of the dark gas belonging to atomic gas. 

7. Conclusions 

We used the Planck-HFI and IRAS data to determine all sky 
maps of the thermal dust temperature and optical depth. The 
temperature map traces the spatial variations of the radiation 
field intensity associated with star formation in the Galaxy. This 
type of map is very important for the detailed analysis of the dust 
properties and their spatial variations. 

We examined the correlation between the dust optical depth 
and gas column density as derived from H i and CO observations. 
These two quantities are linearly correlated below a threshold 
column density of N^^^ < 8.0 x lO^^Hcm"^ corresponding to 
Ay < 0.4 mag. Below this threshold, we observed dust emis- 
sivities following a power-law with ^ 1.8 below A ^ 500 /im 
and flattening at longer wavelengths. Absolute emissivity values 
derived in the FIR are consistent with previous estimates. 

This linear correlation also holds at high column densities 
(N^^^ > 5 X lO^^Hcm"^) corresponding to Ay = 2.5 mag 
where the total column density is dominated by the molecu- 
lar phase for a given choice of the Xco factor. Under the as- 
sumption that the dust emissivity is the same in both phases, 
this leads to an estimate of the average local CO to H2 factor 
of Xco = 2.54 X 1020H2cm-V(Kkms-^). The optical depth in 
the intermediate column density range shows an excess in all 
photometric channels considered in this study. We interpret the 
excess as dust emission associated with dark gas, probably in the 
molecular phase where H2 survives photodissociation, while the 
CO molecule does not. 

In the solar neighbourhood, the derived mass of the dark gas, 
assuming the same dust emissivity as in the H i phase is found to 
correspond to ^ 28%% of the atomic mass and ^ 1 18%% of the 
molecular gas mass. The comparison of this value with the recent 
calculations for dark molecular gas around clouds more mas- 
sive than the ones present in the solar neighbourhood indicates a 
dark gas fraction about three times larger in the solar neighbour- 
hood. The threshold for the onset of the dark gas transition is 
found to be ^ 0.4 "^and appears compatible to, although slightly 
larger than, the thresholds predicted by this model. Finally, we 
stress that the Hi 21 cm line is unlikely to be optically thin and 
to measure all the atomic gas. Therefore, the dark gas detected 
here could well represent a mixture of dark molecular and dark 
atomic gas seen through its dust emission. For an average Hi 
spin temperature of 80 K, the mixture is predicted to be 50% 
atomic and 50% molecular. 

Acknowledgements. A description of the Planck Collaboration and a list 
of its members can be found at http://www.rssd.esa.int/index.php? 
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